Comment on the "Create Dirac Cones in Your Favorite Materials", by Chia-Hui Lin 

and Wei Ku (.arXiv: 1303.4822) 
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Recently a paper by Lin and Ku 1] was posted, where the authors propose a very interesting 
idea to engineer the Dirac points in material which originally did not have such points by inducing 
a CDW state through introduction of impurities like vacancy, substitution, or intercalation. In this 
comment we would like to explain the appearance of the Dirac points in such engineered structures 
by symmetry arguments. 

PACS numbers: 73.22.Pr 



Recently a paper by Lin and Ku [H was posted, where 
the authors propose a very interesting idea to engineer 
the Dirac points in material which originally did not have 
such points. This is achieved by inducing a CDW state 
through introduction of impurities like vacancy, substi- 
tution, or intercalation. In this comment we would like 
to explain the appearance of the Dirac points in such 
engineered structures by symmetry arguments. 

The authors of Ref. [l| actually propose to substi- 
tute the original atoms uniformly across the system in 
such a way, that a lattice with the basis is obtained. We 
will show below that Dirac points are feasible in any lat- 
tice with the basis which possess the space-time inversion 
symmetry The consideration throughout the com- 
ment is in the framework of a tight-binding model with 
the Hilbert space containing one orbital per atom. 

To support our statement let us recall first a well 
known case in which the basis contains two atoms. Hence 
the Hamiltonian at the Bloch wave vector, k = (/ci, fc2), 
is expressed by a 2x2 Hermite matrix, which can be ex- 
panded as 



i?(k) =i;o(k)/ + R(k)-a, 



(1) 



where EQ{\<i) and R(k) = (it!i(k), i?2(k), i?3(k)) are are 
real functions of k, and (cri, (T2) cs) are the Pauli matrices. 
The contact of two bands, described by the Hamiltonian 
([1]), takes place at k, where 



R(k) = 0. 



(2) 



We have two unknowns, ki and ^2, and thus Eq. ([2|) 
is overdetermined. To make the point contact feasible, 
a constraint is required, which reduces the number of 
conditions by one. Such a constraint may be supplied by 
a symmetry. 

Let the Hamiltonian i?(k) possess the space-time in- 
version symmetry, that is 



ff(k) = ITH{k)T-^r^ = CTiiJ*(k)CTi, 



(3) 



where / is the inversion operator and T is the time- 
reversal operator. Substituting the Hamiltonian ([T]) into 
Eq. ^ we obtain = 0, which makes the bands con- 
tact feasible. 



The cases considered by Lin and Ku are more compli- 
cated. In particular, they consider two examples of the 
square lattice; in one case, by substitution of 1/3 of the 
original atoms uniformly across the system, a lattice with 
the basis containing 3 atoms is obtained, in another case, 
by substitution of 1/4 of the original atoms, a lattice with 
the basis containing 4 atoms is obtained (Figs. 3a and 
4a of Ref. [l[ respectively) . 

However, the simple mathematics presented above and 
treating the basis with two atoms, can be generalized for 
the case of any basis containing ns atoms and possessing 
the space-time inversion symmetry Q . The point is that 
instead of direct diagonalization of the Hamiltonian ma- 
trix i?(k), one could renormalize it into smaller dimen- 
sion following the formalism of Brillouin and Wigner. If 
the ns dimensional Hilbert space S is divided into m di- 
mensional subspace Sa and ns — m dimensional subspace 
Sb^ the ns x ns Hamiltonian matrix and its resolvent can 
be written as 



^ ^ / ilAAik) HAB{k) 
^ ' \HBA{k) HBB{k) 



(4) 



(5) 



The resolvent matrix in the subspace Sa can be written 
as 



GAA{k,e)^{e-H''^f{k)) 



(6) 



with the energy dependent effective Hamiltonian pre- 
sented by an m X ni matrix 



H''ff{k,€) = HAA{k) 

+Hab (k) \e - Hbb (k)l ' Hba (k 



(7) 



The full spectrum of the Hamiltonian ^ can be found 
from the energy dependent effective Hamiltonian ([7]). 

Looking at the Figs. 3a and 4a of Ref. [l[ we realize 
that in both cases the basis possess the inversion sym- 
metry, and the inversion is accompanied by the exchange 
of one pair of atoms; other atoms in the basis are not 



2 



moved by the inversion (of course, two atoms separated 
by a lattice vector we treat as the same one) . Taking the 
Hilbert space Sa as that spanned by the orbitals of 
the pair of atoms exchanged by the inversion, we obtain 
Eq. ([3]) in a modified form 



H(k) = ITHCk)T-^I-^ = aiH*(k)ai, 



(8) 



where 




(9) 



for the basis presented on Fig. 3a, and 



^0 1 0^^ 

10 

10 

V 1 / 



(10) 



for the basis presented on Fig. 4a. Thus Haa satisfies 



Eq. (121), Hbb satisfies equation 

ITHBB{'k)T-'r' = H*bb(M), 
and Hab and Hba satisfy equations 

ITHABCk)T-^r' = (TiH*ab{^) 



ITHBA{k)T-^r^ 



(11) 



(12) 
we 



Substituting Eqs. ([2]), (HH) and ^ into Eq. 
recover the two-atoms basis case considered in the be- 
ginning of the comment, that is H'^-f^ {k, e) itself satisfies 
Eq. @. 

Thus we connect the existence of the Dirac points in 
the square lattice examples of Ref . [3] with the fact that 
Dirac points are feasible in any lattice with the basis 
which possess the space-time inversion symmetry. 

Discussion with Chia-Hui Lin is gratefully acknowl- 
edged. 
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